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Abstract 

Numerical investigations around a transformation of Landau's formula 
suggest certain statistical regularities in the distribution of zeros of the 
Riemann zeta function. 

1 Landau's formula 

Denote by Z the set of non-trivial zeros of the Riemann zeta function 
Z = {peC: C(p) = 0, < 9ft(p) < 1} 

and let 

S={-i{p-\):p^Z} 

(so that the Riemann hypothesis is just the assertion that SCM). The standard 
counting function for the zeros is 

N(T) = #(Sn[0,T]) 

and its growth is 

7V(T) = ^log^-^ + 0(logT) (1) 

Z7T ZTT ZTT 

as T — > oo. 

Landau's formula is a statement about the asymptotic behaviour of the 
function 

A -( T ) = E eiat 

tes 

\t\<T 

for general a > 0. Notice that X a (T) is real, regardless of the Riemann hypoth- 
esis, since S = —S and any non-real t € S come in conjugate pairs. 

Landau[5] proved the following theorem. Here A(x) is the von Mangoldt 
function, which is equal to \ogp if x — p k for p prime and k G N + , and 
otherwise. 

Theorem 1 (Landau's formula). For all x > 1, 

^ = 4a(^ + « (2) 

as T — > oo. 
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Thus A a (T) grows linearly with T when a — klogp. By contrast, when 
a ^ klogp its growth is extremely slow: positive and negative values of cos (at) 
for t E S cancel out with great efficiency. 

The ordinates of the zeros of the zeta function around height T share many 
statistical properties with the arguments of the eigenvalues of a random matrix 
from the circular unitary ensemble of order n where n « log ^ . Now it is 
known [I] that the sum of the eigenvalues of such a matrix (all of which lie on 
the unit circle) is, in the limit as n — > oo, distributed as a Gaussian variable with 
fixed variance. In other words, under the random-matrix model the sum of the 
cosines of a sequence of log successive zeros around height T is statistically 
of order 1. But even this goes nowhere near explaining the actual degree of 
cancellation. 



2 Sum over a cycle 

Define 

»?«,h(n) = -log-+ £ (e iat -l). 
a a *■ — ' 

t+heS 

LgJ=» 

We have the following transformation of Landau's formula. 
Theorem 2. If x > 1 and h £ R then 

0<n<M b 



M 



Proof. Write a = logx and T = —M. We use equations [T] and [51 and Stirling's 
formula, to deduce 

£ r)a,h{n) = \e- lah \a(T) - N(T) + i lo g r(^) - J log a + O(logT) 

0<n<M 

T , aT T T . 

+ 7T~ 1°S 7. 7J 7T~ 1°§ a + O(logT) 
Z7T Z7T Z7T Z7T 

= -^A(x)x- ( -i+^ + O(logT) 

Z7T 

and the result follows. □ 

Numerical investigations, using data that Odlyzko has made available on 
his website[3], suggest that the distribution of T] a ^(n) around its mean value is 
asymptotically independent of n, as well as of h. 

Conjecture 1. For every a > there is a B or el function f a : C — > M satisfying 

fa = 1 
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such that for every h £ K and every Borel set ICC, 

#{0 < n < M : %ih (n) + Afc!l e -(|+*)a e X} y 
,i im ' U = / 

M— >oo M J x 

It is natural to consider the closely related continuous function H a : R + — ► C 
defined by 

Ha (r) = -log-!-- £ (l + e - fa ). 
a Z7ra * — ' 

ti|<7T 

Here t — aS denotes {t — at :t £ S}. It is easily checked that for all a and h, 

Va , h (n) = H a ((2n + 1)tt + ofc) + 0(0 (3) 

as n — > oo, so that the asymptotic distribution of rj a ^h{n) is determined by that 
of H a (r) for r = 7r + ah modulo 2n. 

As t increases continuously, H a (r) switches between near-circular orbits, as 
illustrated in Figure [TJ At 'time' r the orbit is centred on the point 

-l°g^- - v a {r) 
a Zira 

where 

v a (T) = #((T-aS)n[-n,n}). 

These centres lie on the real axis at unit spacings, and they drift slowly 
in the positive direction. However, it seems that the statistical distribution of 
H a (r) is asymptotically independent of t. 

Conjecture 2. For every a > there is a Borel function F a : C — > R satisfying 

F n = 1 



such that for every Borel set ICC, 



v -x T 1 



X 



where A denotes Lebesgue measure. 



Relation [3] implies that the function F a must be related to the function f a 
in Conjecture [T] by 

F a {z) = ±- f a (z-^e-i°+ i(> )d9. 
2vr J a 

In particular, the two functions are the same when a ^ fclogp. 

Figures [21 [3] and 2] illustrate the (empirical) density functions F a for a = 1 , 
a = | and a = log 2 respectively. Figure [5] illustrates /i og 2- 

These plots are based on calculations from the first 100000 zeros. The dis- 
tribution seems to be independent of height, and plots based on zeros around 
the (10 12 ) th , (10 21 ) st and (10 22 ) nd zero look similar. But even if the conjectures 
should in the end prove false, the patterns evident in these plots do demand an 
explanation. 
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Figure 1: Locus of H x (t) for 50000 < r < 50600. 
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Figure 2: Values of H\(t) in the complex plane, for 50000 values of t equally 
spaced between it and 74920 — n. 
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Figure 4: Values of H\ og 2(r) in the complex plane, for 50000 values of t equally 
spaced between n and 74920 log 2 — n. 
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